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ABSTRACT 

A search has been made for neutrinos from the hep reaction in the Sun and from the diffuse 
supernova neutrino background (DSNB) using data collected during the first operational phase 
of the Sudbury Neutrino Observatory, with an exposure of 0.65 kilotonne-years. For the hep 
neutrino search, two events are observed in the effective electron energy range of 14.3 MeV 
< T e ff < 20 MeV where 3.1 background events are expected. After accounting for neutrino 
oscillations, an upper limit of 2.3 x 10 4 cm _2 s _1 at the 90% confidence level is inferred on the 
integral total flux of hep neutrinos. For DSNB neutrinos, no events are observed in the effective 
electron energy range of 21 MeV < T e ff < 35 MeV and, consequently, an upper limit on the v e 
component of the DSNB flux in the neutrino energy range of 22.9 MeV < E v < 36.9 MeV of 
70 cm" 2 s _1 is inferred at the 90% confidence level. This is an improvement by a factor of 6.5 on 
the previous best upper limit on the hep neutrino flux and by two orders of magnitude on the 
previous upper limit on the v e component of the DSNB flux. 

Subject headings: neutrinos, Sun: general, supernovae: general 
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1. Introduction 

The Sudbury Neutrino Observatory (SNO) is 
a real-time, heavy water Cherenkov detector lo- 
cated in the Inco, Ltd. Creighton nickel mine near 
Sudbury, Ontario, Canada at a depth of 6010 m 
water equivalent (Boger et al. (2000)). SNO de- 
tects electrons and neutrons from, respectively, 
charged-current (CC) and neutral-current (NC) 
interactions of neutrinos on deuterons, as well as 
neutrino-electron elastic scattering (ES) interac- 
tions, in one kilotonne of D2O contained in a 
12 m diameter acrylic vessel (AV). These inter- 
actions are observed via Cherenkov light detected 
by 9456 photomultiplicr tubes (PMTs) mounted 
on a 17.8 m diameter support structure. By com- 
paring the observed rates of these interactions, 
SNO has demonstrated that a substantial fraction 
of the 8 B electron neutrinos produced in the Sun 
transform into other active neutrino flavors (Ah- 
mad et al. (2001, 2002a,b); Ahmed et al. (2004); 
Aharmim et al. (2005)). These results are con- 
sistent with the predictions of neutrino oscilla- 
tions (Maki, Nakagawa & Sakata (1962); Gribov & 
Pontccorvo (1969); Wolfenstein (1978); Mikheyev 
& Smirnov (1985)). 

The Sun generates energy by nuclear fusion; 
protons combine to form helium in reactions that 
release neutrinos. The primary solar fusion pro- 
cess is a series of reactions known as the pp chain. 
Five reactions in the pp chain produce neutri- 
nos; the highest energy neutrinos are those from 
the hep reaction: 3 Hc + p — > 4 He + c+ + v e . 
The endpoint of the hep neutrino spectrum is 
18.77 MeV and lies above that of the 8 B spec- 
trum, which is approximately 15 MeV. The flux 
of hep neutrinos (e.g. Bahcall & Krastev (1998), 
Bachall & Pinsonneault (2004)) is currently pre- 
dicted to be (7.97 ± 1.24) x 10 3 enr^- 1 (Bahcall, 
Serenelli & Basu (2005)), which is small compared 
to the fluxes from the other neutrino-producing 
reactions in the pp chain, including the 8 B flux 
which has been measured to be (4.95 ± 0.42) x 
10 6 cm^s" 1 (Aharmim et al. (2005)). The dom- 
inant contribution to the uncertainty in the hep 
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neutrino flux prediction is 15.1% from the calcula- 
tion of the nuclear matrix elements (Park et al. 
(2003)). The previous best upper limit on the 
hep neutrino flux is 7.3 x 10 4 cm _2 s _1 at the 
90% confidence level (CL), based on measurements 
with the Super-Kamiokande detector (Hosaka et 
al. (2005)). After accounting for neutrino oscil- 
lations, this limit can be interpreted as an up- 
per bound on the total flux of hep neutrinos of 
1.5 x 10 5 cm s . Currently, only one reaction 
( 8 B) from the pp chain has been uniquely observed 
and measured experimentally. An observation of 
hep neutrinos would give further confirmation of 
the pp chain as the primary solar energy genera- 
tion mechanism and would allow further tests of 
the solar model. 

Neutrinos produced in core-collapse supernovae 
also contribute to the energy region above the 8 B 
endpoint. The current generation of neutrino de- 
tectors can detect the transient signal from a su- 
pernova in the Milky Way, but the expected signal 
from a supernova in a more distant galaxy is fewer 
than one event. Neutrinos from all cxtragalactic 
supernovae since the beginning of the formation of 
stars in the universe constitute the diffuse super- 
nova neutrino background (DSNB), which may be 
detectable. Model predictions range from 0.19 to 
1.49 cm _2 s _1 for the v e component of the DSNB 
flux in the neutrino energy range 22.9 MeV < 
E v < 36.9 MeV (Beacom & Strigari (2006), Ando 
& Sato (2003)). The best upper limit on the 
77 e component of the DSNB flux is 1.2 cm _2 s _1 
at the 90% CL for E„ > 19.3 MeV, based on 
measurements with the Super-Kamiokande detec- 
tor (Malek et al. (2003)). While an indirect limit 
on the v e component of the DSNB flux can be in- 
ferred from this (Lunardini (2006)), the previous 
best direct upper limit is 6.8 x 10 3 cm _2 s _1 for 
neutrino energies 25 MeV < E v < 50 MeV, based 
on measurements with the Mont Blanc liquid scin- 
tillator detector (Aglietta et al. (1992)). 

A search for hep and DSNB neutrinos has been 
performed by counting the numbers of events in 
predefined energy intervals (signal boxes) above 
the 8 B endpoint. The most sensitive signal boxes 
for this analysis were selected by evaluating the 
predicted signal and background levels before ex- 
amining the data. Given the predicted signal and 
background levels in the signal boxes, limits on 
the flux of hep and DSNB neutrinos are set us- 



ing a modified Fcldman-Cousins technique. The 
following sections describe the data set, detector 
response, determination of the backgrounds, anal- 
ysis procedures and the limits obtained for the hep 
and DSNB neutrino fluxes. 

2. The Data Set 

The data included in these analyses were col- 
lected during the initial phase of SNO operation, 
during which the detector contained pure D 2 0. 
The data were collected from 1999 November 2 
until 2001 May 28 and comprise 306.4 live days 
corresponding to an exposure of 0.65 kilotonne- 
years (Ahmad et al. (2002a)). 

Since results from this phase were last pub- 
lished, numerous improvements have been made 
to the analysis tools, many of which were used in 
the analysis of data from phase two (Aharmim et 
al. (2005)), for which two tonnes of salt were dis- 
solved in the heavy water. Further improvements 
were applied in this analysis, the most significant 
of which was improved estimation of the effective 
electron kinetic energies (T c ff ) of the events, based 
on the optical paths to each operational PMT. 
Other enhancements include improved handling 
of false hits due to cross-talk between electronics 
channels and an improved accounting of working 
PMTs using both neutrino and calibration data 
to track bad channels. However, the vertex recon- 
struction algorithm was the same as that used in 
previous phase one analyses, in which events were 
reconstructed under the assumption that they are 
due to single electrons. This is more suited for the 
reconstruction of hep and DSNB events than the 
algorithm used in phase two. After the application 
of the new analysis tools, events inside the kinetic 
energy window of 12 MeV < T e gf < 35 MeV were 
not examined until the hep and DSNB signal boxes 
had been selected. 

In addition to the event selection discussed 
in Aharmim et al. (2005), which includes a se- 
lection that removes Michel electrons with visible 
precursors, selection criteria were applied to re- 
move backgrounds from atmospheric neutrino in- 
teractions. As the hep and DSNB signals are ex- 
pected to be single electron events, these back- 
grounds can be reduced significantly by remov- 
ing events which correlate in time with neutrons, 
electrons or 7-rays. Consequently, any candidate 
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event that appeared within 250 ms of another 
with T e ff > 4 MeV and a reconstructed vertex 
inside the AV was removed. In addition, two 
Kolmogorov-Smirnov (KS) tests were applied: one 
to test the azimuthal symmetry of the PMT hits 
about the reconstructed event direction; the other 
to test the compatibility of the angular distribu- 
tion of PMT hits with that expected from a single 
electron. In the signal boxes, the selections on 
PMT hit isotropy and the prompt light fraction 
were further tightened with respect to previous 
SNO analyses (Aharmim et al. (2005)). This was 
possible in this analysis due to the higher energies 
of the candidate events. The combined event selec- 
tion reduced the expected number of atmospheric 
neutrino events in the hep signal box by a factor of 
29 and by a factor of 77 in the DSNB signal box. 
The signal acceptance of the combined event selec- 
tion is (96.6 ±0.7)% for hep and (94.0 ±1.5)% for 
DSNB events, measured using calibration source 
data and simulation. 

3. Detector Response 

To understand the signals and backgrounds 
in this analysis, it is important to measure the 
energy response and uncertainties in the signal 
boxes. The energy response can be parameter- 
ized by a Gaussian of resolution er T = — 0.154 + 
0.390-/7T + 0.0336T e , where T e is the true kinetic 
energy of the electron. In SNO analyses, Monte 
Carlo simulation is used to estimate the response 
of the detector to different particles. The prop- 
agation of electrons, positrons and 7-rays is car- 
ried out using EGS4 (Nelson, Hirayama & Rogers 
(1985)). The uncertainties in the energy scale 
and resolution of the SNO detector have typi- 
cally been measured using 6.13 MeV 7-rays from 
a 16 N source (Dragowsky et al. (2001)). At the 
higher energies more characteristic of this analy- 
sis, Michel electrons from muon decays and a pT 
( 3 H(p,7) 4 Hc) source (Poon et al. (2000)), which 
produces 19.8 MeV 7-rays, were used to comple- 
ment the 16 N measurements. Using simple event 
selection criteria, including one based on the time 
between events, 135 Michel electrons were iden- 
tified in the data. Potential deviations in energy 
scale and energy resolution between data and sim- 
ulations were assumed to be linear functions of 
energy. These functions were fit with a maximum 
likelihood technique using data from 16 N and pT 



sources as further constraints. The results were 
used to refine the energy scale and resolution es- 
timates and to measure their uncertainties at the 
analysis thresholds. An energy scale uncertainty 
of 0.96% and a resolution uncertainty of 3.8% were 
estimated at the hep threshold of 14.3 MeV. At the 
DSNB threshold of 21 MeV, an energy scale un- 
certainty of 1.06% and a resolution uncertainty of 
6.0% were estimated. Correlations between these 
quantities were included in the final analysis. Ad- 
ditional non-Gaussian tails to the resolution func- 
tion were also considered, but were found to be 
insignificant. Data and Monte Carlo distributions 
of T e ff for 16 N and pT calibration events and for 
Michel electrons are shown in Figure 1. 

Event vertex and direction reconstruction were 
unchanged from the analysis in Ahmad et al. 
(2002a). The position resolution at 15 MeV is 
(12.0 ± 2.5) cm and the angular resolution is 
(20.6 ± 0.4)°. These were measured using a com- 
bination of 16 N source data and simulation. The 
same fiducial volume, defined by events recon- 
structed within a distance of 550 cm from the cen- 
ter of the detector, was selected. The uncertainty 
on the expected number of events within the fidu- 
cial volume due to vertex accuracy was 2.9%. 

4. Backgrounds 

Three distinct classes of background are consid- 
ered: 8 B neutrino interactions, atmospheric neu- 
trino interactions and instrumental backgrounds. 
Figure 2 shows the simulated energy spectra of 
the signals and backgrounds, normalized to their 
expected rates. 

Electrons from 8 B neutrino interactions are the 
dominant (97%) background for the hep analysis 
but are a negligible background for DSNB. These 
events can reconstruct into the hep signal box due 
to the finite energy resolution of the detector. The 
magnitude of the 8 B background depends on the 
details of the detector response and is very sensi- 
tive to the energy scale and resolution at thresh- 
old. In the CC interaction, by which SNO predom- 
inantly detects the 8 B and hep neutrinos, there is 
a strong correlation between neutrino and electron 
energy. This, in addition to a cross section that 
rises with the square of the energy rather than 
linearly, provides a clearer distinction between the 
two neutrino spectra in the region of the 8 B end- 
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Fig. 1. — The effective electron kinetic energy 
spectra from data and Monte Carlo for (a) events 
from the 16 N source, (b) events from the pT source 
and (c) Michel electrons. The data are shown 
in the energy regions free of source related back- 
grounds. 
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Fig. 2. — The simulated effective electron kinetic 
energy spectra of the signals and backgrounds of 
interest in the hep and DSNB analyses. Also 
shown are the data in the range 6 MeV < T or f < 
12 MeV which are used to normalize the 8 B elec- 
tron and neutron distributions. The atmospheric 
neutrino background is made up of a number of 
different signals: neutrons at low energies; 7-rays 
from nuclear de-excitations at intermediate ener- 
gies; and Michel electrons and CC interactions of 
atmospheric u e and u e on deuterons at higher en- 
ergies. The DSNB model in this figure is the 
T = 6 MeV model from Beacom & Strigari (2006). 
The third class of background, instrumental back- 
grounds, which is not shown in this figure is neg- 
ligible. 
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point than is possible with the ES interaction. 

The 8 B background also depends on the de- 
tails of the shape of the detected electron spec- 
trum. The 8 B neutrino spectrum from Winter 
et al. (2003, 2006) was assumed along with its 
quoted uncertainties. Neutrino oscillations were 
taken into account by correcting and combining 
the electron spectra from CC and ES interactions 
using the energy-dependent v e survival probabil- 
ity from the joint solar neutrino and KamLAND 
(Araki et al. (2005)) oscillation analysis presented 
in Aharmim et al. (2005). Additional spectral ad- 
justments were included to account for CC inter- 
actions on 18 O, radiative corrections to the CC 
deuteron cross section (Nakamura et al. (2002); 
Kurylov et al. (2002)) and the acceptance of the 
event selection. The expectation for the hep signal 
is constructed in the same way using the neutrino 
spectrum from Bahcall & Ulrich (1988) with cor- 
rections from Bahcall (1997). 

After the determination of the 8 B signal shape, 
its normalization was determined using data at 
lower energies, where the hep signal is insignifi- 
cant. In an energy window of 6 MeV < T ff < 
12 MeV, 2006 events were observed. To determine 
normalizations, these data were fit using a maxi- 
mum likelihood technique with probability density 
functions (PDFs) for the 8 B electrons (CC and 
ES signals) and neutrons (NC signal and back- 
ground). The distributions used in this fit were 
functions of event energy and direction with re- 
spect to the Sun (cos 9q). The results of this fit 
were then used to estimate the 8 B contribution 
inside the signal boxes. 

Atmospheric neutrino interactions produce a 
second class of background events. They are the 
dominant background in the DSNB signal box and 
come from several sources: 

• Electrons from low energy (E v < 100 MeV) 
charged-current v e and F e interactions; 

• Michel electron events, where the precursor 
muons (and pions) are below the Cherenkov 
threshold and do not trigger the detector; 

• 15.1 MeV 7-rays from de-excitation of an ex- 
cited state of 12 C created via a nuclear cas- 
cade from neutrino interactions on 16 O; 

• Mis-identified non-electron events. 



For low energy atmospheric v e and V e , the flux 
prediction from Battistoni et al. (2005) is used, 
which has an uncertainty of 25%. Only charged 
current interactions on deuterons, with cross sec- 
tions from Nakamura et al. (2002) and Kurylov et 
al. (2002) are considered; the contributions from 
other interaction types are not significant. The in- 
teractions of these neutrinos constitute 14% of the 
DSNB background, but are insignificant in the hep 
signal box. 

Events from the latter three sources are associ- 
ated with atmospheric neutrinos of higher energy 
(E u > 100 MeV). Monte Carlo simulations were 
used to generate atmospheric neutrino interactions 
in the SNO detector with statistics equivalent to 
500 times the expected number of events. For this 
purpose, the package NUANCE (Casper (2002)) 
was used, with the Bartol04 flux prediction for 
Sudbury (Barr et al. (2004)). The flux uncertainty 
in the neutrino energy range that contributes to 
the background is 10%. The events generated by 
NUANCE were then propagated and fully simu- 
lated in the SNO Monte Carlo, from which back- 
ground predictions were obtained after application 
of the event selection. 

To assess uncertainties these events were di- 
vided into three categories. The first category, 

quasi-elastic (QE) CC events, is the primary 
source of untagged Michel electrons and origi- 
nates from neutrinos in the energy range 150 to 
250 MeV. The uncertainty on the cross section in 
this energy region is 25% (Barish et al. (1977)). 
These Michel electrons comprise 80% of the DSNB 
background. For the second category, 15.1 MeV 7- 
ray events, there is no data in the literature on pro- 
duction rates and, thus, a 100% uncertainty was 
assigned to the production rate predicted by NU- 
ANCE, which uses the calculation of Ejiri (1993). 
These 7-rays constitute half of the atmospheric 
background in the hep analysis, but due to the 
magnitude of the 8 B background constitute only 
1.5% of the total hep background. The final cat- 
egory comprises QE NC events and interactions 
that produce pions, to which a cross section un- 
certainty of 30% is assigned (Ahrens et al. (1987)). 

There is an additional uncertainty applicable to 
the latter two categories of atmospheric neutrino 
interactions. A comparison of events from data 
and the simulation has shown that the simulation 
underestimates the production of correlated neu- 
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trons. It is unclear if this is due to errors in the 
prediction of primary neutron production or in the 
transport of hadrons in the simulation. However, 
there is good agreement between data and Monte 
Carlo for correlated electron events. Events in the 
simulation are re-weighted in such a way that the 
average neutron multiplicity is changed to better 
match the data. This results in a change to the 
background rejection rate in the simulation due to 
time correlated neutrons. This correction results 
in an additional uncertainty of 7% in the rate of 
atmospheric background events inside the signal 
box that are not due to QE CC interactions. After 
application of the event selection to remove events 
with correlated neutrons, electrons and gamma- 
rays the atmospheric background in these analyses 
is reduced by a factor of two. 

To verify the predictions for the atmospheric 
neutrino background, data outside the signal box 
in the energy range 35 MeV < T e g < 55 MeV 
were examined. This energy range was selected 
to be most sensitive to the main component of 
the atmospheric neutrino background: the Michel 
electrons. In this energy range, 0.28 Michel elec- 
trons and 0.05 electrons from low energy charged- 
current atmospheric neutrino interactions are ex- 
pected. One event was observed, consistent with 
the predictions of the simulation. Inside this en- 
ergy range, the effect of the event selection on 
events correlated with neutrons, electrons or 7- 
rays was also examined. Two such events were 
observed, each consistent with being an otherwise 
untagged Michel electron preceded by a 7-ray from 
the de-excitation of the nucleus participating in 
the primary neutrino interaction. No events were 
observed that are correlated with neutron or elec- 
tron events. These results are also consistent with 
the predictions of the simulation. 

The final class of backgrounds is associated with 
instrumental effects such as electronic pickup or 
static discharge from the PMTs. For these events, 
an upper limit of 0.002 events is set in an energy 
range of 6 MeV < T c ff < 35 MeV using the tech- 
nique described in Aharmim et al. (2005). This 
analysis is not sensitive to the isotropic acrylic 
vessel background (IAVB) (see Aharmim et al. 
(2005)). To predict the number of IAVB events 
that pass the signal event selection, the 13 IAVB 
events clearly identified in the data by simple cri- 
teria arc scaled via Monte Carlo simulation. A 



limit of 7 x 10 4 IAVB events in the energy range 
of 14 MeV < T off < 35 MeV is inferred at the 90% 
CL. 

5. Analysis and results 

The analysis was designed to construct confi- 
dence intervals on the neutrino fluxes using a mod- 
ified Feldman-Cousins approach (Cousins & Feld- 
man (1998); Conrad et al. (2003); Hill (2003)). 
Limits were also calculated using a Bayesian ap- 
proach (Eidelman et al. (2004)); very similar re- 
sults are obtained for the two techniques. To de- 
termine confidence limits the probability p(N\S) 
of observing N events, given a signal flux S, 
is calculated taking statistical fluctuations and 
all known systematic uncertainties into account. 
A Monte Carlo technique is used to integrate 
over the estimated distributions of the system- 
atic uncertainties, including known correlations, 
by sampling ensembles of shifted parameter val- 
ues and propagating their effect on the PDFs and 
extracted signal and background normalizations. 
The major uncertainties included in this proce- 
dure are shown in Table 1. 

The hep and DSNB analyses are very similar, 
except that the definitions of signal and back- 
ground arc modified. In the DSNB analysis, the 
hep distribution is scaled using the standard solar 
model prediction, including its uncertainty, and 
added to the background estimate. 

For the hep analysis the signal box was cho- 
sen to optimize the sensitivity based on Monte 
Carlo simulations. The sensitivity was defined as 
the mean value, from an ensemble of Monte Carlo 
experiments, of the 90% CL upper limit for the 
hep flux, integrated over all energies using the hep 
neutrino spectrum and accounting for neutrino os- 
cillations as discussed in Section 4, assuming the 
standard solar model. Figure 3a shows the pre- 
dicted numbers of signal and background events 
with their la uncertainties as the lower threshold 
of the hep signal box is changed, and Figure 3b 
shows the sensitivity of the analysis as a function 
of the signal box threshold. There is a region be- 
tween 12.5 MeV and 14.3 MeV where the sensitiv- 
ity is nearly flat. Within this range the choice of 
the best signal box is a compromise between the 
signal to background ratio and signal acceptance. 
The energy range 14.3 MeV < T eff < 20 MeV was 
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Table 1: Major uncertainties included in the anal- 
yses. 



Source of Uncertainty 


Magnitude of effect 


Energy scale 




l e fT = 14. 6 MeV 


U.yo/o 


J- cff — ^i ivie v 


1 .UD /o 


Energy resolution 




l e ff = 14. o MeV 


Q 00/ 

0.070 


i off — -^i ^ie V 


r no/ 

O.U/o 


Vertex accuracy 


z.y /o 


Vertex resolution 


2.5 cm 


A n en l mi' T'ficn 1 1 1" i nn 


2% 


8 B v e spectrum 


Taken from 1 


tan 2 (6>i 2 ), Ato 2 2 


Contours from 2 


V a tm flux 




E v > 100 MeV 


10% 


E v < 100 MeV 


25% 


Cross sections 




CC deuteron 


1.2% 


Vatm CC QE 


25% 


v atm other 


30% 


15.1 MeV 7-rays 


100% 


v atm n-multiplicity 


7% 



References. — (1) Winter et al. (2006); (2) Aharmim et 
al. (2005) 
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Fig. 3. — Panels (a) and (b) show the expected 
number of events and the sensitivity of the analy- 
sis respectively as the energy threshold of the hep 
signal box is varied. The hep flux in panel (a) 
is normalized to the solar model prediction. The 
upper limit of the signal box is fixed at 20 MeV. 



selected. In this range, the variations of predicted 
signal and background levels due to systematic un- 
certainties are strongly correlated, as can be seen 
in Figure 4. In this signal box, 3.13 ± 0.60 back- 
ground events and 0.99 ±0.09 signal events are ex- 
pected. The contributions to the signal and back- 
ground uncertainties from the dominant sources of 
systematic uncertainties are shown in Table 2. 

Two events are observed in the hep signal box. 
After accounting for the effect of neutrino oscilla- 
tions, this results in an upper limit on the integral 
total hep neutrino flux of 2.3 x 10 4 cm _2 s _1 at 
the 90% CL. This is 2.9 times the prediction of 
the standard solar model. Using a Bayesian tech- 
nique rather than the modified Feldman- Cousins 
approach, a limit of 2.9 x 10 4 cm~ 2 s _1 is found 
at the 90% CL. The spectrum of events in the re- 
gion of the signal box is shown in Figure 5. The 
shape agrees with the background prediction at 
the 77.8% CL based on Monte Carlo simulations 
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Fig. 4. — The joint probability distribution of sig- 
nal and background events in the hep signal box 
assuming the standard solar model flux of hep neu- 
trinos. 
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Table 2: The effect of systematics on hep signal 
and background 



Systematic Source 


ABkg(%) 


ASig(%) 


TT in ormr cr*£i 1 
_1_j11L.I Jiy stale 


1 % 7 


7 


Energy resolution 


9.7 


0.7 


Vertex accuracy 


0.3 


2.9 


8 B v e spectrum 


0.8 


0.0 


v atm flux 


0.3 


0.0 


Am 2 2 


0.6 


0.5 


tan 2 (0i2 ) 


0.7 


3.2 


Cross sections 






CC deuteron 


0.0 


1.1 


Vatm CC QE 


0.3 


0.0 


15.1 MeV 7-rays 


0.8 


0.0 


Low energy fit statistics 


3.1 


0.0 


Combined Width 


19.1 


9.0 



Note. — This tabic shows the one standard deviation con- 
tributions to the width of the signal (ASig) and background 
(ABkg) probability distributions in the hep signal box. The 
combined widths are greater than the quadrature sums of 
the systematics due to correlations and non-linearities. 



using a KS statistic. 

This result is model-dependent, as neutrino os- 
cillations were assumed in the predictions of sig- 
nal and background. A model-independent search, 
in which a limit is placed on the v e component 
of the integral hep flux, was also carried out. 
This search is independent of any model of neu- 
trino flavor change. A further event selection of 
cos 00 < 0.8 was applied to remove ES events 
(which are directed away from the Sun) and leave 
events due only to CC (v e ) interactions. The en- 
ergy threshold was selected so that the mean back- 
ground expectation was less than 0.25 events, and 
flux limits were set conservatively by assuming 
there is no background and quoting only the up- 
per bound. With these criteria, a signal box of 
16 MeV < T cff < 20 MeV was selected. Without 
neutrino oscillations 0.66 ± 0.08 hep events are ex- 
pected in this box. As no events are observed a 
limit on the v e component of the integral hep neu- 
trino flux of 3.1 x 10 4 cm^s" 1 at the 90% CL is 
inferred. 

Signal box selection for the DSNB search fol- 
lows the approach of the hep neutrino search. A 
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Fig. 5. — The distribution of events in the region 
of the 8 B endpoint. There are two events in the 
hep signal box 14.3 MeV < T cff < 20 MeV. Also 
shown are the estimated number of background 
events, including the systematic uncertainty, and 
the standard solar model prediction for the hep 
signal. 



signal box of 21 MeV < T e ff < 35 MeV was se- 
lected in this case. This choice was bounded at 
the high end by the prior choice of the hidden en- 
ergy interval and at the low end by a desire to 
minimize any contribution from hep neutrinos. In 
this signal box, 0.18 ± 0.04 background events are 
expected. The expected number of signal events 
depends upon the assumed DSNB model, but is 
uncorrelated with the background prediction as 
different systematic uncertainties are dominant for 
signal and background. The estimated values of 
the systematic uncertainties are shown in Table 3. 
No events are observed in the signal box, result- 
ing in an upper limit of 2.3 events due to DSNB 
neutrinos at the 90% CL. 

To obtain a DSNB flux limit, a spectral model 
for the DSNB neutrinos is required. In this paper, 
three models from Beacom & Strigari (2006) and 
two models from Ando & Sato (2003) for differen- 
tial flux predictions have been used. Table 4 shows 
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Table 3: The effect of systematics on DSNB signal 
and background 









Energy scale 


3.4 


0.9 


Energy resolution 


1.7 


0.4 


Vertex accuracy 


2.9 


2.9 


Event selection 


0.9 


1.7 


v fluxes 






Vatm E v > 100 MeV 


8.6 


0.0 


Vatm E v < 100 MeV 


3.5 


0.0 


hep 


0.2 


0.0 


Cross sections 






Vatm CC QE 


19.2 


0.0 


v atm other 


2.6 


0.0 


CC dcuteron 


0.2 


1.2 


v a tm n-multiplicity 


0.7 


0.0 


Combined Width 


24.7 


4.3 



Note. — This tabic shows the one standard deviation con- 
tributions to the width of the signal (ASig) and background 
(ABkg) distributions in the DSNB signal boxes. For the 
DSNB signal the T = 6 MeV model from Beacom & Strigari 
(2006) is used. As in the prediction for the signal and back- 
ground in the hep signal box, combined widths are greater 
than the quadrature sums of the systematics. 

the integral flux predictions for these models and 
the 90% CL upper limits inferred from data. 

Using these results, a limit can also be derived 
on the DSNB v e flux for neutrinos that produce 
electrons with kinetic energies inside the DSNB 
signal box. Although the integral flux upper lim- 
its are significantly different for these models, since 
their spectral shapes are similar in the signal box, 
the resulting upper limits for the neutrinos in this 
region vary little (see Table 4) . Taking the median 
result a limit on the DSNB v e flux of 70 cm _2 s _1 
at the 90% CL for 22.9 MeV < E v < 36.9 MeV is 
inferred. These limits and the background predic- 
tion are in good agreement with those predicted 
by Beacom & Strigari (2006), after accounting for 
the difference in exposure between their prediction 
and the data used in this search. 

6. Conclusions 

Data from the first operational phase of SNO, 
with an exposure of 0.65 kilotonne-years, have 



been used to search for neutrinos from the hep re- 
action in the Sun. No evidence for these neutrinos 
was observed, and an upper limit on the integral 
total flux of hep neutrinos of 2.3 x 10 4 cm~ 2 s -1 has 
been derived at the 90% CL. This measurement 
improves the previous best limit on the hep neu- 
trino flux, measured with the Super-Kamiokandc 
detector, by a factor of 6.5, but is not inconsis- 
tent with the standard solar model. A model- 
independent limit on the integral hep v e flux, with 
no assumptions about neutrino oscillations, is set 
at 3.1 x 10 4 cm _2 s _1 . A search for the v e compo- 
nent of the diffuse supernova neutrino background 
has also been made using SNO data. Again, no ev- 
idence for these neutrinos was found, and an upper 
limit at 90% CL on the v e component of the DSNB 
flux of 70 cm-V 1 for 22.9 MeV < E v < 36.9 MeV 
is inferred. This is an improvement of two orders 
of magnitude on the previous v e limit (Aglietta et 
al. (1992)). The exposure of the final SNO data 
set for these analyses combined across all phases of 
the experiment, is expected to be approximately 
four times that of the data used in this analysis. 
A future search for hep and DSNB fluxes using 
these data will be carried out, which is expected 
to further improve upon the limits presented in 
this paper. 
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Table 4: DSNB flux predictions and limits 







Integral Flux 


Flux 22.9 MeV < E v < 36.9 MeV 






(cm 2 s v ) 




(cm 2 s - 1 ) 




Model 




Prediction 


Upper Limit 


Prediction 


Upper Limit 


B&S : 


T = 4 MeV 


21.1 


1.1 x 10 4 


0.19 


93 


B&S : 


T = 6 MeV 


14.1 


1.5 x 10 3 


0.66 


72 


B&S : 


T = 8 MeV 


10.5 


6.0 x 10 2 


1.08 


61 


A&S : 


NOR-L 


28.5 


1.3 x 10 3 


1.49 


69 


A&S : 


NOR-S-INV 


34.9 


2.3 x 10 3 


1.06 


70 



Note. — This table shows the 90% CL upper limits on the u e component of the DSNB flux and model predictions for different 
models from Beacom & Strigari (2006) (B&S) and Ando & Sato (2003) (A&S). 
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